function [dat]=r_seawifs
% ---- Read SeaWiFS parameter
% (chlorophyll concentrations)

% ---- Create the Structure ...        
dat = struct('dummy',1);

% ---- NOCS
%[dat2]=load_gadget([str2num(year) (str2num(year)+0.9999999999)],[-104 52 -21 82],'chlo_sw8D_v10a.nc');
%dat.londata=dat2.lon;
%dat.latdata=dat2.lat;
%dat.timedata=dat2.time;
%dat.data=dat2.chlo;
% From decimal year to MATLAB Serial number date
%dat.timedata=datenum(fix(dat.timedata),1,1)+((dat.timedata-fix(dat.timedata)).*(datenum(fix(dat.timedata),12,31)-datenum(fix(dat.timedata),1,1)+1));

%dat.data=permute(dat.data,[3 2 1]);

% ---- LEGOS
% sdir='/calcul/pc-charria/OBSERVATIONS/SATELLITE/OC/';
% liste=dir(sdir);
% nn=char(liste.name);
% tagdate=nn(3:end,28:35); % Last day of the week
% dat.timedata=datenum([tagdate,repmat('T000000',[size(tagdate,1) 1])],'yyyymmddTHHMMSS');
% %dat.data(1:380,1:160,1:size(tagdate,1))=NaN;
% h=waitbar(0,'Read and reshape seawifs fields ...');
% for ific=1:size(tagdate,1)
%     waitbar(ific/size(tagdate,1))        
%     temp=load([sdir,nn(ific+2,:)]);
%     dat.londata=unique(temp(:,1));
%     dat.latdata=unique(temp(:,2));
%     dat.data(1,1,ific)=temp(1,3);
%     rla=temp(1,2);
%     ilo=1;
%     ila=1;
%     for il=2:size(temp,1)
%         if (rla==temp(il,2))
%             ilo=ilo+1;            
%         else
%             ilo=1;
%             ila=ila+1;
%             rla=temp(il,2);            
%         end
%         dat.data(ilo,ila,ific)=temp(il,3);
%     end
%         
% end
% close(h)
% 
% dat.data(dat.data > 900)=NaN;


% ---- NOCS - from Toulouse
sdir='/calcul/pc-charria/OBSERVATIONS/SATELLITE/OC/';
load([sdir,'Sea_Nemo.mat'])
dat.londata=lon_chl;
dat.latdata=lat_chl;
dat.timedata=time_chl;
dat.data=chl;
% From decimal year to MATLAB Serial number date
dat.timedata=datenum(fix(dat.timedata),1,1)+((dat.timedata-fix(dat.timedata)).*(datenum(fix(dat.timedata),12,31)-datenum(fix(dat.timedata),1,1)+1));

dat.data=permute(dat.data,[3 2 1]);

% ---- Clean the Structure ...
dat=rmfield(dat,'dummy');    
 